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Abstract 

I consider the Diophantine approximation problem of sup-norm simultane¬ 
ous rational approximation with common denominator of a pair of irrational 
numbers, and compute explicitly some pairs with large approximation con¬ 
stant. One of these pairs is the most badly approximable pair yet computed. 

The theory of approximation of a single irrational number by rationals is well 
known, and for our purposes the relevant facts may be summarized as follows. 
We measure the goodness of approximation of the rational number p/q to a by 
c(a,p,q ) = q\qa—p\. For each irrational a (without loss of generality, we may 
assume 0 < a < 1) we know by Dirichlet’s thereom that there are infinitely many 
rationals p/q such that \a—p/q\ < l/q 2 , or c(a,p, q) < 1. It is therefore of interest 
to ask how small one may make yin c(a,p. q) < ybefore this property fails to hold. 
The approximation constant of a is thus defined as c(a) = liming —KX> c(a,p,q). 
Here, of course, for each q we choose the p whch minimizes c(a,p,q). Numbers 
a with a large c(tt) are hard to approximate by rationals. The one-dimensional 
Diophantine approximation constant, defined as c { — sup a , ]R c(a'), has the value 
I/a/ 5, attained at a — (V5 — 1)/2. 

Otherwise expressed, this means that c l is the unique number such that for 
each £ > 0, the inequality c(a,p, q) < c : +£ has infinitely many rational solutions 
p/q for all a, whereas there is at least one a such that c(a,p, q) < c l —£ has only 
finitely many rational solutions. 
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These results completely solve the problem of rational approximation in one 
dimension, but by contrast the situation in two or more dimensions is much more 
complex and in fact the value of the analogous constant c„ for 2 is unknown 

B0- 

We wish to simultaneously approximate a pair of irrationals by a pair of ra- 
tionals with common denominator and to measure the closeness of approximation 
by the maximum error in the two components, so we make the definitions: for 

P = (P\■,P2) e ^ Z ,a = (oc v a 2 ) G M 2 , let 

c(a,p,q) — <7 max ( \qCX l —p l \ 2 , \qa 2 —p 1 \ 2 ) 

and 

c(a ) = liminf { c(a,p,q), p G Z 2 , q G Z }. 

q —>00 

The two-dimensional (sup-norm) simultaneous Diophantine approximation con¬ 
stant is then 

c 2 = sup c(a). 

aeR 2 

Despite much work over the last few decades the value of c 2 

is unknown, though folk-lore suggests that its value is 2/7. Adams [4] has shown 
that this is the correct value if we restrict the pair (oCj, cx 2 ) to cubic number fields, 
but his result does not give us a constructive procedure to identify pairs with large 
c(a). 

Here, however, I use a theorem of Cusick together with high-precision nu¬ 
merical computation to explicitly compute examples of such pairs. These have 
potential applications to numerical simulation studies of dynamical systems on 
the 2-torus, where (dCj, cc 2 ) represent the winding number of periodic orbits. 

Cusick’s construction makes use of the cubic number field Q(0), where 0 — 
2cos(2;r/7), of smallest positive discriminant, namely 49. For details on cubic 
fields and their integral bases, I refer to |Q]. 

The theorem of Cusick [fj states that for any integral basis (1, ce,/3} of Q (0), 
we have c*< 2/7, where c* is the infimum of those c such that 

\x+ay+Pz\ max(y 2 ,z 2 ) < c 

(with y and z not both zero) has infinitely many solutions in integers x,y,z. Addi¬ 
tionally, for any e > 0 there is an integral basis {1, a, /3} such that 

2/7—c*(a,/3) < £ 
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iff 


1: The continued fraction of 0 has patterns [... ,n x , 1, l,n 9 ,...] with n { ,n 2 ar¬ 
bitrarily large; or, 

2: The continued fraction of 0 has patterns [...,n v 2, n 2 ,.. .] with n l ,n 2 arbi¬ 
trarily large. 

It is not known whether either of the last two conditions are satisfied. Note that 
this theorem relates to the dual problem to simultaneous Diophantine approxi¬ 
mation, namely approximation to zero by linear forms. Hence, it is not imme¬ 
diately apparent that the upper bound of 2/7 that it gives for c* it defines is rel¬ 
evant to the problem of determining c 2 . However, from another paper by Cu- 
sick ([]9p, Corollary 1 on page 187), we have that for the particular field Q(0), 
c*(a,/3) = c(a.p) for all integral bases. Also, by a theorem of Davenport [|T(j|]. 
we have supc*(a,/3) = supc(a,/3), where the sups are over all irrational pairs, 
not necessarily in a cubic field. 

Thus, if the above patterns in the continued fraction of 0 do in fact exist, 
Cusick’s theorem gives us a way of finding explicit pairs (which together with 1 
form an integral basis of Q( 0 )) with a value of c close to 2/7. Even if n x , n 2 do not 
become arbitrarily large, just the presence of some large values gives us potential 
candidates for very badly approximable pairs. 

From results in [|6|], it follows that for an integral basis of the form {l,pO + 
q0 2 1 r0+s0 2 },( p r ^) e PSL( 2, Z), where —q/p and —s/r are rational approximants 
to 0 obtained by truncating the continued fraction at the points where condition 1 
or 2 is satisfied, c* is explicitly given by 

c* = 1/max{ |A+5+C|, \A—B+C\, \C—B 2 / (4A)|, \A—B 2 /(AC) \ } 

= 1/max{ \A+B+C\, |A-5+C|,49/|4A|,49/|4C| }, 

where 
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with 


a = (0 2 2 -0 2 )(0 2 -0 1 2 ) 

b = (e 2 2 -e 2 )(e l -e)+(02-e)(ef-« 2 ) 

c = (e-e 2 )(e t -e) 

0 = 2cos(2;r/7) 

0j = 2cos(4;r/7) 

0 2 = 2cos(6;r/7). 

With this background, I can now state the main result of this paper: I have 
exactly computed over 2 million partial quotients of the continued fraction of 0 
(directly from the defining cubic x 3 +x 2 — 2x— 1), and the required patterns do 
indeed occur, though very infrequently. The largest values of c*, with the corre¬ 
sponding fractional parts of a — pO+qO 2 and /3 = r0+s6 2 occur at: 

(A) positions 57-60: [...,60,1,1,50,...], c* w 0.2851877 

a w 0.4563286858107963651609830446124431560745665647128596153008802 
f3 w 0.4781573193903170892895817415258772866671562381178937772663665 

(B) positions 2927-2930: [...,22,1,1,22,...], c* w0.2853154 

a w 0.1554011929520066325796747316744656830061413509133865038820677 
j8 w 0.6003679362632065361061389158735863615694126556922931077332356 

(C) positions 3626-3629: [...,272,1,1,215,...], c* « 0.2855726 

aw 0.6530646111210617321254054547968773238346090082060701183776580 
/3w 0.9410463762107594592302548739412493098027738320829952592216557 

(D) positions 33877-33880: [...,81,1,1,78,...], c* w 0.2856261 

a w 0.9319638477108390366188499907354642637920661848031694636081724 
J3 w 0.7032571495109702868148790086182835032528572663181375225766851 

(E) positions 215987-215990: [...,124,1,1,129,...], c* w 0.2856678 
aw 0.4375520476578757564544576313180510209212270982522655674846137 
jSw 0.5646614639128419094417646922292433724548272488193131214134926 
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(F) positions 957740-957743 [...,460,1,1,415,...], c* « 0.28568046 
a « 0.6134980317071692745070006892224661159462079954445253478668675 
J8« 0.9411544329571988683307282702558980820407618535628393885417987 

(G) positions 1650050-1650053: [...,648,1,1,666,...], c* ^0.2857082 
« -0.4848739572889332951989678247806190621159456336657613155291560 
P « 0.5404925035004667478257428539575752367424111926723566428410541 


These calculations involve extremely large integer and floating-point numbers; 
in case (G) the absolute values of the integers p,g,r, s are of the order 2 3xl ° 6 , 
and the calculation of c* requires floating-point operations of about twice this 
precision. In fact, these examples all come from cases of Cusick’s first condition, 
and c* is given by 49/ |4A | or 49/ |4C|. Of course, the approximate decimal values 
for a,P given above are insufficient to represent the true values, but these may be 
reconstructed if required from the continued fraction of 0. 

An independent verification of these results may be obtained by giving the val¬ 
ues a,/3 as input to a simultaneous Diophantine approximation algorithm. Such 
an algorithm finds all best simultaneous approximants up to a given denominator. 
For the computation of sup-norm best approximants, an algorithm has been given 
by Furtwangler [ ]TT| , [Hj ]. Figure [I] shows the behaviour of the Furtwangler algo¬ 
rithm applied the pairs (A) and (C) above. The approximation constant estimated 
from the minimum c after ignoring the initial transient is about 0.2856, verifying 
the more precise value of c* above. But the chief point to be noted is the extremely 
long initial transient. Until a sufficient large denominator q is reached, these pairs 
would in fact appear to be not badly approximable. 

I have thus exhibited some explicit pairs which are very badly approximable 
by rationals. I believe that the value 0.2857082 above is the largest explicitly com¬ 
puted lower bound for the two-dimensional simultaneous Diophantine approxima¬ 
tion constant c-,. 

The question remains open as to whether there are pairs (necessarily unrelated 
to the field Q(0)) with approximation constant larger than 2/7. 
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Figure 1: c(a.(3,q) vs. log 10 (g) at best approximants for two integral bases 

(l,ce,j3) of the field Q(0). Above: case (A), below: case (D). 
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